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1.  Introduction  and  summary 

We  consider  stationary,  additive,  interval  functions  X(a).  These  are  vector 
valued  stochastic  processes  having  real  intervals  a  =  (a,  /?]  as  domain,  having 
finite  dimensional  distributions  invariant  under  time  translation  and  satisfying 

(1.1)  X(AjUA2)  =  X(a  i )  +  X(a2), 

for  disjoint  intervals  Aj.  a2.  Such  processes  are  considered  in  some  detail  in 
Bochner  [5].  Setting 

(1.2)  X(/)  =  X(0, /], 

—  oc  <  /  <  oo,  and  in  the  reverse  direction  setting 

(1.3)  X(a,/J]  =  X(0)  -  X(a), 

we  see  that  a  consideration  of  stationary  interval  functions  is  equivalent  with  a 
consideration  of  processes  X(t),  —  oo  <  t  <  oo,  having  stationary  increments. 
These  last  are  discussed  in  Yaglom  [24]  for  example.  Important  examples  of 
processes  of  the  type  under  consideration  are  provided  by  the  point  processes. 
Here  the  components  of  X(a)  give  the  number  of  events  of  various  sorts  that 
occur  in  the  interval  a.  A  variety  of  properties  and  applications  of  point  pro¬ 
cesses  may  be  found  in  Cox  and  Lewis  [11],  Bartlett  [4],  and  Srinivasan  [21]. 

The  paper  is  divided  into  various  sections.  In  Section  2  we  introduce  a  key 
assumption  for  the  processes;  specifically  we  require  that  all  the  moments  of 
X(a)  exist  and  have  particular  integral  representations.  We  are  then  able  to 
define 


(1-4)  fai, •••,«<< (^1  ’  >  ^k)> 

—  oo  <  kj  <  oo,  «!,*•*,  ak  =  1,  •  •  •  ,  r,  the  cumulant  spectra  of  order  k  of  the 
r  vector  valued  X(a).  These  turn  out  to  be  generalizations  of  the  cumulant  spectra 
of  order  &  of  a  continuous  time  series  discussed  in  Brillinger  and  Rosenblatt 
[9].  We  then  present  a  spectral  representation  for  X(a).  This  representation  was 
introduced  in  Kolmogorov  [17]  for  real  valued  processes  with  stationary  incre¬ 
ments.  It  takes  the  form 


(1.5) 


X(0,  t] 


*00 

exp  [iXt]  —  1 

J  —  00 

iX 

^Zx(/l), 
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with  Zx{X),  —  oo  <  X  <  go,  a  vector  valued  process.  The  process  Zx(X)  relates 
to  the  cumulant  spectrum  (1.4)  through  the  expression 

(1.6)  cum  {dZai(Xx  ),•••,  dZau{Xk)} 

=  ..aMw- ,di.k 

with  3(X)  denoting  the  Dirac  delta  function. 

In  Section  3  of  the  paper  we  indicate  how  the  theory  developed  applies  to 
integrated  continuous  time  series,  to  point  processes,  and  to  processes  that  are 
hybrids  of  these  last  two.  In  the  case  of  point  processes  we  relate  the  cumulant 
spectra  to  important  parameters  that  have  been  introduced  by  Bahba  [1]  and  by 
Kuznetsov  and  Stratonovich  [18]. 

Section  4  of  the  paper  discusses  various  asymptotic  properties  of  the  statistic 

(1.7)  d*/^)  =  f  exp  { —  iXt}  dX(t) 

Jo 

based  on  an  observed  stretch  of  an  X(a)  process.  It  will  be  seen  to  behave  in  a 
similar  manner  to  the  finite  Fourier  transform  of  a  stretch  of  a  continuous 
stationary  series.  It  follows  that  the  estimates  of  the  various  cumulant  spectra 
of  X(a)  may  be  formed  in  the  manner  of  Brillinger  and  Rosenblatt  [9]  and  that 
the  properties  developed  in  that  paper,  such  as  asymptotic  normality,  will  con¬ 
tinue  to  hold.  A  selection  of  results  that  therefore  become  available  is  provided. 
In  particular  results  relating  to  the  linear  time  invariant  regression  of  one 
stationary  interval  function  on  another  are  given.  Because  point  processes  are 
particular  cases  of  the  processes  under  consideration  it  follows  that  an  asymptotic 
theory  for  the  spectral  estimates  of  order  two  of  point  processes  has  now  been 
provided. 

In  Section  5  we  apply  the  previously  mentioned  asymptotic  results  to  develop 
estimates  of  the  parameters  suggested  by  Bahba  and  by  Kuznetsov  and  Stratono¬ 
vich  for  point  processes.  In  Section  6  we  consider  the  problem  of  the  estimation 
of  the  second  order  spectra  of  a  continuous  time  series  when  its  values  are  avail¬ 
able  only  for  random  times  that  are  the  occasions  of  events  of  an  independent 
point  process. 

Section  7  discusses  briefly  some  practical  implications  and  extensions  of  the 
results  of  the  paper.  The  proofs  of  the  various  lemmas  and  theorems  of  the  paper 
are  given  in  Section  8. 

I  would  like  to  thank  P.  A.  W.  Lewis  for  a  variety  of  helpful  comments  on  the 
point  process  sections  of  the  paper. 


2.  Random  interval  functions 

Let  A  denote  the  collection  of  finite  intervals  of  the  form  a  =  (a.  /?].  We 
consider  r  vector  valued  stochastic  processes  X(a),  a  g  A  with  the  additivity 


property 

(2.1) 
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X(A)  =  X( A ! )  +  X(A2), 

for  A.  Aj .  a2  €  A  with  a  =  Aj  u  a2,  At  n  a2  —  0.  Such  a  process  will  be  called 
an  r  vector  valued  additive  stochastic  interval  function.  As  one  example  we  mention 

(2.2)  X( a )  =  J  Y(/)  dt. 

where  Y(£),  —  go  <  t  <  oo,  is  a  continuous  r  vector  valued  time  series.  As  a 
second  example  we  consider  X(a)  =  N(a)  where  N(a)  is  an  r  vector  valued  point 
process  with  N0(a)  giving  the  number  of  events  of  the  ath  component  of  the 
process  that  occur  in  the  interval  a.  As  a  final  example  we  mention 

(2.3)  X(a)  =J  Y(t)dN(t). 

where  Y(t),  —  go  <  /  <  oc,  is  an  r  vector  valued  continuous  time  series  and 
At(a).  A  G  A,  is  a  point  process.  If  tj,  •  •  •  ,  t„  denote  the  times  of  events  of  the 
process  At(a)  in  the  interval  a.  then  X(a)  equals 

(2.4)  Y(t,  )  +  •  •  •  +  Y(t„), 
in  this  case. 

In  connection  with  the  process  X(a),  a  g  A,  we  set  down 
Assumption  2.1.  The  process  X(a),  a  g  A,  is  an  r  rector  valued  stochastic 
interval  function  possessing  moments  of  all  orders  such  that  for  Aj.  •  •  •  .  Ak  e  A; 
ai>  ‘  '  ‘  >  ak  =  1,  •  •  •  ,  r\  k  =  1,  2,  •  •  •  , 

(2.5)  E{X.^)  ■  ■  ■  Xn(*k)}  =  f  ■  f 

J Ai  J Ak 

for  some  function  Mau. . .  ak(tl,  •  •  •  ,  tk),  —  oo  <  tj  <  oc.  of  bounded  variation  in 
finite  intervals. 

In  the  case  that  X(a)  satisfies  this  assumption  and  f>a{t)  is  bounded  and  con¬ 
tinuous  for  t  in  some  interval  of  A  and  0  outside  the  interval,  we  may  define 
stochastic  integrals  of  the  form 

(2.6)  J  dXJt) 

as  the  limit  in  mean  (of  any  order  v  >  0)  of  approximating  Riemann  sums 

(2.7)  t  X„(A j). 

j=  1 

where  tj  g  A^andAi  u  •  •  •  u  a„  is  a  partition  of  the  support  of  </)a(0, a  =  1,  *  •  •  ,  r. 
(See  Cramer  and  Leadbetter  [12],  p.  86  for  the  case  v  =  2.)  These  integrals  have 
the  property 
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(2.8)  fijj  4v(f.)  <«„,((,)  J  <j>n(lk)dXn(tk)  | 

=  J  •  •  •  J  <t>c  (<l)  ‘  •  •  <l>ak(h)  dMal  ...  ak(t1,  •  •  •  ,  tk), 

for  ax,  •  •  •  ,  ak  =  1,  •  •  •  ,  r. 

For  a  =  (a,  /?]  in  A  we  denote  the  translated  interval  (oc  +  t,  +  t]  by  a  +  t 
for  —  co  <  t  <  oo.  We  will  now  say  that  an  r  vector  valued  additive  stochastic 
interval  function  is  stationary  if  the  joint  distributions  of  all  finite  collections  of 
variates 

(2.9)  X(At  4-  t),  •  •  •  ,  X(Afc  +  t), 

a !,*•*,  Aft  6  A,  —  go  <  t  <  oo,  and  k  =  1,  2,  •  ■  •  ,  do  not  depend  on  t.  In  this 
connection  we  have 

Lemma  2.1.  If  X(a),  a  e  A.  is  a  stationary  r  vector  valued  additive  interval 
function  satisfying  Assumption  2. 1,  then  for  al,  •  •  •  ,  ak  =  1,  •  •  •  ,  r\k  =  1,  2,  •  •  •  , 

(2.10)  EiX^A,)  ■  ■  ■  XnI.Ak)} 

=  f  f  dMau  •  •  •  ,ak(h  ~  hi  '  ' '  ’  h-  i  —  h)  dtk, 

J  A  j  %)  Afc 

for  some  function  Mf  ...  ak(u1;  ••  •  ,uk_1),  —  oc  <  Uj  <  oc ,  of  bounded  variation 
infinite  intervals. 

In  the  case  k  =  1 .  the  lemma  indicates  that 

(2.11)  EXa(  a)  =  c;|a|, 

for  some  constant  C'a .  a  =  1,  •  •  •  ,  r  with  |  a  |  denoting  the  length  of  the  interval  a. 
It  follows  from  this  lemma  that  one  can  write 

(2.12)  cum  {XJAj  ),•••,  XajAk)} 

=  r  r  dcai'...,ak(tl  —  tk,  ■  •  • ,  tk_x  —  tk) dtk, 

«/  A  i  %)  Aj< 

for  a  function  C'ai  ...  ak(ux,  •  •  •  .  uk_  1),  of  bounded  variation  in  finite  intervals. 
In  differential  notation  we  may  write  this  last  as 

(2.13)  cum  {dXai(ul  +<).*•*.  dXak_fuk~x  +  /),  dX0k(t)} 

=  dC'au...  ak(u i,  •  •  •  ,  uk_1)dt. 

Taking  note  of  Lemma  2.1  and  (2.12)  we  set  down  the  key  assumption  of  our 
work.  It  is 

Assumption  2.2.  The  process  X(a),  a  g  A  .is  a  stationary  r  vector  valued  addi¬ 
tive  interval  function  satisfying  Assumption  2.1  and  such  that  C'a  1  ...  ( u  j,  •  •  • ,  uk_  , ) 

of  (2.12)  satisfies 

(2-14)  /_•••/ \u}\  d\C'ai  ...  ak(ul,  ■  •  •  ,  uk  —  j )  |  <  oo, 

for  j  —  ,  k  :  a1?  •  •  •  ,  ak  =  1,  •  •  •  ,  r  ;  k  =  2,  3,  •  •  •  . 
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This  assumption  has  the  nature  of  a  mixing  condition  on  the  increments  of 
X(/),  that  is,  increments  that  are  well  separated  in  time  are  only  weakly  dependent. 
In  view  of  condition  (2.13)  we  can  define  the  Fourier  transforms 


(2-15)  ll5*--,Ak) 

=  (271  )~k+l  j  ■■■Jexpj-*£  t-JU]\dC'au---.ak(Ul>  ’  ‘  '  >  Uk-l)< 


for  —  oo  <  Xx,  •  •  •  ,  Xk  <  oo,  where  we  understand  Xx  +  •  •  •  +  Xk  =  0  in  the 
definition.  For  completeness  we  set 

(2.16)  fjX)  =  C„ 

where  C’a  was  defined  in  (2.11),  a  =  1,  •  •  •  ,  r.  The  transform  fau.. .  ,ak(^i ,  ’  '  ’  ,  Xk) 
is  called  a  cumulant  spectrum  of  order  k  of  the  process  X(a),  a  g  A.  We  will 
sometimes  find  it  convenient  to  adopt  the  unsymmetric  notation 

(2.17)  .,a„  •  •  • ,  xk., )  =  „k(A,,  •  ■  ■ .  xk). 

The  second  order  cumulant  spectra,  ffb( A).  —  oo  <  X  <  oo,  are  of  particular 
importance.  It  is  convenient  to  collect  them  together  into  the  r  x  r  spectral 
density  matrix 

(2-18)  l'x,x(^)  =  [/a,b(^)]‘ 

We  also  collect  the  first  order  spectra  together  into  the  r  vector 

(2.19)  rx  =  [/;]. 

There  is  an  intimate  connection  between  stationary  interval  functions  and 
stationary  series.  Suppose  that,  X(a),  a  g  A,  satisfies  Assumption  2.2.  and  has 
cumulant  spectra 


(2.20)  ...(A,,  ••,4). 

Suppose  the  real  valued,  (f)a{t ),  —  oo  <  t  <  go,  satisfies 


(2.21)  J \t\\(pa(t)\ dt  <  go, 

fora  =  1,  ■  •  •  ,  r.  Then  the  r  vector  valued  times  series,  Y(/),  —  oo  <  t  <  oc.with 
components 


(2.22) 


Ya{t)  =  J  f)a{t  -  u)  dXa(u), 


a  =  1,  •  •  •  ,  r  may  be  seen  to  be  stationary  and  such  that 
(2.23)  cum  {Yai(t  +  /,),  •  •  •  ,  Yak_ft  +  tk_,),  Yak(t)} 

=  J  •  •  •  J  4>atf  1  -  «l)  ’  ’  ‘  <l>ak-Mk-  1  -  Uk- l)  4>aMk) 

C,,.  .,«k(«C  *••»«*-!)  duk ■ 
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Taking  a  Fourier  transform,  we  see  that  the  cumulant  spectra  of  Y (t),  in  the 
sense  of  Brillinger  and  Rosenblatt  [9],  are  given  by 


(2.24)  .  U 
where 

(2.25)  =  J  exp  {- iXt }  <j>a{t)dt, 
for  a  =  1,  •  •  •  ,  r. 

A  variety  of  authors  (including  Kolmogorov  [17],  Doob  [14]  p.  551,  Ito 
[16],  Yaglom  [24],  [25]  p.  86,  Bochner  [5]  p.  159)  have  given  spectral  repre¬ 
sentations  for  stationary  interval  functions  (or  processes  with  stationary  incre¬ 
ments).  In  this  connection  we  mention 

Theorem  2.1.  Let  the  process  X(a),  a  g  A,  satisfy  Assumption  2.2.  Let 


(2.26) 


Z«/'(A)  =  (2k)~' 


f T  1  —  exp  {  —  iXt] 

J-r  —  it 


dX(t), 


for  —  oo  <  X  <  oo.  Then  there  exists,  Zx(X),  —  oo  <  X  <  oo,  such  that  Z(P{X) 
tends  to  Zx(X)  in  mean  order  v,  for  any  v  >  0.  Zx(X)  satisfies 


(2.27)  cum  {Zai(Xt),  •  •  •  ,  Zait(Xk)} 


•  •  •  ,  a*)  dal  ■  •  •  dak, 


for  al}  •  '  ,ak  =  1,  •  •  •  ,  r\  k  =  1,  2,  •  •  •  .  Also 

(2.28)  X(a)  =  i:[i  exp  [iXt]  d/J  dZx{X), 
with  probability  one. 

In  differential  notation  particular  cases  of  (2.27)  include: 

(2.29)  EdZx(X)  =  d{X)fxdX\ 

(2.30)  Cov  {dZx(X),  dZx(p)}  =  5(X  -  y)tx,x(X)  dX  dp; 

(2.31)  cum  {dZafXx),  •  •  •  ,  dZak(Xk)} 

—  d(Xx  4-  •  •  •  +  Xk)fau...  ak(X1,  ■  -  •  •  ,  Xk.l)  dXx 

Also  if  we  set  X(t)  =  X(0,  /],  then  (2.28)  takes  the  form 

exp  {iXt}  — 


dX  t 


(2.32) 


x(«)-f[“p{u}~1]  ”,(*)■ 


for  —  oo  <  t  <  oo . 
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The  representation  (2.28)  is  useful  for  displaying  the  effect  of  linear  time  in¬ 
variant  operations  on  the  process,  X(a),  a  g  A.  Suppose  a(A),  a  g  A,  is  an  s  x  r 
matrix  valued  interval  function  of  bounded  variation  satisfying 

(2.33)  J  \t\  d|a(J)|  <  oo. 

Set 

(2.34)  A(A)  =  Jexp  {  —  iXt)  da,(t), 

for  —  oo  <  X  <  oo .  The  s  vector  valued  interval  function 

(2.35)  Y(A)=a*x(A) 

a(A  —  u)  dX(u), 

a  e  A,  may  be  seen  to  satisfy  Assumption  2.2.  Also  the  process  Zy(A),  —  oc  < 
X  <  oo,  of  its  spectral  representation  may  be  seen  to  satisfy, 

(2.36)  dZY(X)  =  A(A)  dZx(X). 

We  may  infer  from  this  last  that  the  spectral  density  matrices  of  X(a)  and  Y(a) 
are  related  by 

(2.37)  fy,y(A)  =  •A.(A)fxix(^)A(A)  . 

This  last  relation  has  the  identical  form  with  that  giving  the  effect  of  linear  time 
invariant  operations  on  the  spectral  density  matrices  of  time  series. 

We  conclude  this  section  by  remarking  that  the  function  M'au...ak{u  i,  *  •  • ,  uk_  x) 
of  (2.10)  may  be  determined  as 

(2.38)  M'au...tak(u1,  •  •  •  ,  ttkii) 

T -1  fT  xai(t,t  +  Ml]  •**Xak_1(«,«  +  Mk_i]dXflk(<) 

Jo 


3.  Some  examples 

Example  3.1.  Suppose  that  Y(£),  —  oo  <  t  <  oo,  is  an  r  vector  valued 
stationary  time  series  possessing  moments  of  all  orders.  If 

(3.1)  c'au...tak{ux,  •  •  •  ,  uk_l)  =  cum  {Yai{t  +  ut),  •  •  •  ,  Yak^(t  +  Mk_j),  Yak(t)}, 

satisfies 


(3.2) 


J  •  "  J  |«Jc;if...tflk(M1,  •  •  •  ,  wk_i)|  dux  ■  ■  •  duk_i  <  00. 
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the  cumulant  spectra  of  the  series,  Y (t),  — oo  <  t  <  oo,  are  given  by 
(3-3)  '  '  '  >  A-k) 


understanding  +  •  •  •  +  Xk  =  0  (see  Brillinger  and  Rosenblatt  [9]).  Also  the 
Cramer  representation  of  Y (t)  is  given  by 

(3.4)  Y  (t)  =  J  exp  {iXt}  dZY(X), 


where 


(3.5) 


ZY(X)  =  l.i.m.  (2n)  1 
r- oo 


exp  {  —  iXt]  —  1 
—  it 


Y (/)  dt. 


Suppose  we  construct  the  interval  process 


(3.6)  X(a)  =  J  Y(<)  dt, 

then  we  quickly  see  that  this  process  satisfies  Assumption  2.2  with 

ru  i  ptih-i 

(3.7)  C';i,...,Bk(a1,--*,ttk_1)  =  J  ---J  c'au...tak(ui,  •  •  •  ,  «*_!) 

dux  ■  •  •  duk_1. 

The  cumulant  spectra  of  the  interval  process,  X(a),  a  e  A,  are  therefore  the  same 
as  the  cumulant  spectra  of  the  time  series  Y  (t),  —cc<t<co. 

A  comparison  of  expressions  (3.5)  and  (2.26)  indicates  that  ZX(A)  of  the 
spectral  representation  of  X(a)  is  equivalent  with  Zy(A)  of  the  Cramer  represent¬ 
ation  of  Y (t). 

In  a  later  section  we  will  see  that  our  proposed  empirical  analysis  of  the  process, 
X(a),  a  g  A,  reduces  to  the  usual  empirical  analysis  of  the  continuous  series 
Y (t),  —oo  <  t  <  oo. 

Example  3.2.  Consider  an  r  vector  valued  point  process,  N(a),  a  g  A.  Here 
Na{ a)  represents  the  number  of  events  of  the  ath  sort  that  occur  in  the  interval  A. 
If  we  let  la  denote  a  vector  with  1  as  its  ath  component  and  0  elsewhere,  then 
we  may  set  down 

Assumption  3.1.  The  point  process,  N(a),  a  g  A,  possesses  moments  of  all 
orders  and  is  such  that  if  a1?  •  •  •  ,  Ak  are  disjoint  intervals  with  | ax |,  •  •  •  | Afc|  ^ 
<5  <  oo, 

(3.8)  />{N(i1)  =  n„  •  •  • ,  N(i»)  =  n,}  <  i:J|41|1"'1  ■  ■  • 


for  some  finite  Kd  and  for  nl9  *  •  •  ,  nk  having  nonnegative  integral  coordinates. 
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Also  if  tjE  Aj,  for  such  a1s  •  •  •  ,  a  k,  there  is  a  f unction  pau...ak(tx  •  •  •  ,  tk), bounded 
infinite  intervals,  such  that 

(3.9)  lim  |41|-1-"|4l|-1i>{N(i,)  =  l«,.-",N(Al)  =  lak} 

|A,|-0 

Pai,  •  •  • ,  »  ’  ^k)» 

uniformly  in  tx,  •  •  •  ,  tk. 

The  functions  pai,---,ak(h>  '  '  ’  •  h)  have  been  called  product  density  functions, 
see  Sriniyasan  [21].  The  function  pa(t),  —  oo  <  t  <  co,  is  called  the  density  of 
events  of  the  ath  sort  at  time  t,  a  =  1 ,  ■  ■  •  ,  r. 

We  note  that  the  process  satisfies 

(3.10)  i’ll  N.W>  l}  =  0(N2), 

and  is  therefore  orderly  (events  tend  not  to  happen  simultaneously).  Also  we 
have 

(3.11)  P{N.(  a)  =  1}  =j).(«)|a|  +  o(|a|), 
and 

(3.12)  P{Na( a)  =  0}  =  1  -  p.(<)|a|  +  o(|a|). 

In  the  theorem  below  we  let  <5{;r}  denote  the  Kronecker  delta  S{x}  =  1  if 
x  =  0  and  <5{a:}  =  0  otherwise.  We  let  XA(x)  denote  the  indicator  function 
Za(t)  =  1  if  t  €  a,  xa(t)  =  0  otherwise.  We  have 

Theorem  3.1.  Let  the  r  vector  valued  point  process  N(a),  a  £  A,  satisfy 
Assumption  3.1.  Then 

(313) 

js{m..(a,)  •  ■  •  Naj Ak)>  =  t  i  cn  «{*i  -  “j}]  •  cn  <5!^  -  <*j}] 

1  at,  •  •  •  ,a(  =  1  jevi  je\( 

f-  ••  fen  wti)]  • "  cn  xajMIpiu-.**  (*i,  ••  •  ,  t,)^  ••  • 

J  J  jevl  jev( 

with  the  sum  extending  over  all  partitions  (vl5  •  •  •  ,  \e)  of  the  set  (1,  •  •  •  ,  k). 

We  see  that  the  moments  of  N(a),  a  £  A,  have  the  integral  representation 
required  in  Assumption  2.1.  Particular  cases  of  this  theorem  include 

(3.14)  ENa{  a)  =  jpa(x)dx, 

(3.15)  J&{2V-(a1)JV-(a2)}  =  f  f  pa,a(*i,  i2)dx1  dx2  +  f  pa(x)dx, 

JAi  c/  A 2  */AioA2 

(3.16)  ^{^.(A1)^t(A2)}  =  f  f  Pa,b^i^2)dxldx2 

J  Ai  J  A 2 


if  a  b, 
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and 

(3.17)  E{NM)k}  =  I  ■  ■  f  ,.(t i.  ' ' ' .  V)*i  ' ' '  */, 

^  =1  Ja  Ja 

where  «9c’^k)  denotes  a  Stirling  number  of  the  second  kind.  Expression  (3.17) 
was  set  down  by  Ramakrishnan  [20]  and  Kuznetsov  and  Stratonovieh  [18]. 

Kuznetsov  and  Stratonovieh  [18]  remarked  that  it  might  prove  more  useful 
to  consider  the  cumulant  functions 

(3.18)  Qai,- •  ,ak(h’  ‘  ‘  ‘  i  h) 

=  (-!)'■  V  -  1)'-Pajijevt(h’j  e  Vi)  •  •  -  Paj;jeve{tpj  e  **)> 

where  the  summation  extends  over  all  partitions  (vt ,  •  •  •  ,  v,)  of  (1 ,  •  •  •  ,  k).  These 
functions  have  the  property  of  tending  to  0  as  |/,  —  tj\  ->  oo  in  the  case  that  the 
increments  of  N(<)  are  tending  to  become  independent  as  they  separate  in  time. 
Particular  cases  of  the  functions  include 

(3-19)  qa{t)  =  pa(t), 

(3.20)  qa,b(t i,  h)  =  Pa.bVu  h)  ~  Pa(h)Pb(h)- 

The  inverse  relation  to  (3.18)  is 


(3.21)  Pai,-  -  ,ak(h> 


h) 


k 

J',  Q.ay,  je  vflj  >  j  ^  V1 ) 


We  have 

Theorem  3.2.  Let  the  r  vector  valued  point  process  N(a),  a  g  A,  satisfy 
Assumption  3.1.  Then 


(3.22) 

cum  {Nai( Aj),  •  •  •  ,  Nak{ Ak)}  =  £  I 

£  =  1  at,  •  ■  • ,  tXf  =  1 

f •  •  •  fen  xa,(ti)]  •  •  •  [n  xa,(^)] 

J  J  jfzv  i  jevf 


cn  ^(ai  -  aji]  •  *  tn  -  ajU 

jev  i  jeve 

‘  '  '  >  xt)  •  •  •  dxe 


where  the  summation  extends  over  all  partitions  (vl5  •  •  •  ,  v^)  of  the  set  (1,  •  •  •  ,  k). 

The  relation  (3.22)  has  the  same  form  as  the  relation  (3.13).  As  particular 
cases  we  mention 


(3.23) 


ENa{ a)  =  qa( t)  dx 


dx, 


(3.24)  Cov  {Na(Ai),  Na(A2)}  =  f  f  Qa,a(zu  x2)dxl  dx2  +  f  qa(z)dx, 

J  Ai  JA  2  JAinA2 


(3.25) 


Cov  {Na(^i),  Arb(A2)}  —  1,  L 


qa,b(z i»  z 2 )  dx i  dx2, 
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and 

(3.26)  cumt  {Na( a)}  =  £  f  •  •  •  f  qa,...,a{ t1?  •  •  •  ,  r,)  dij  •  •  ■  dt,. 

^=1  Ja  Ja 

This  last  expression  was  given  by  Kuznetsov  and  Stratonovich  [18].  We  remark 
that  in  differential  notation,  (3.22)  has  the  form 

(3.27)  cum  {dNai(tt),  ■  ■  ■  ,dNak{tk)}  =  qau...tak(t1,  •  *  •  ,  tk)  dt1  •  •  •  dtk, 

if  the  tj  are  distinct.  As  a  further  implication  of  the  theorem  we  have 

Corollary  3.1.  Under  the  conditions  of  the  theorem  and  if  (j)a(t)  is  continuous 


with  finite  support ,  a  =  1,  •  •  •  , 

k ,  then 

(3.28)  cum  jj  $x(t)  dNai(t), 

•  •  ■  ,  J  4>k(t)  dA7flk(o| 

II 

« 

II 

l - 1 

» 

i  -  «j}]  •  •  •  [fl  3 &  ~  f  ' '  fCn 

jevf  J  J  jev  i 

•••[n  •  •  •  ,  T,)  dxx  "-dxe 

jevf 

where  the  summation  is  over  all  partitions  (Vj .  •  •  •  .  v^)  of  (1 ,  •  •  •  ,  k). 

If  the  point  process  N(a),  a  g  A.  is  stationary,  then 

(3.29)  Pai,- ■  ■  ,ak(t  +  tl,  • 

•  •  •  *  +  tk)  =  Pau...,ak{t1.  •  •  •  •  <*), 

and 

(3.30)  qau--,ak{t  +  <1,  • 

•  •  +  h)  =  .  4)- 

for  all  real  t,  t{,  •  •  •  ,  tk.  In  this  case  we  set 

(3.31)  rau...<ak(u  i,  •  •  •  ,  u*_j)  =  qau..,ak{u  i,  •  .uk_1.  0). 

The  parameter  ra  is  called  the  mean  intensity  of  the  process  A7a(A),  a  e  A:  ra  a(u) 
is  called  the  covariance  density  of  the  process  Na{  a),  A  G  A;  and  rab(u).  fora  =/=  h , 
is  called  the  cross  covariance  density  of  the  component  Na(  a)  with  the  component 
Nb(A). 


We  now  set  down 

Assumption  3.2.  N(a),  a  g  A ,  is  an  r  vector  valued  stationary  point  process 

satisfying  Assumption  3.1  and  such  that 

(3.32)  j...j 

1)|  dux"  -  duk-i  <  go 

for  ax,  •  •  •  ,  ak  =  1,  •  •  •  ,  r  ;  k  =  2,  3,  •  •  ■  . 

If  the  process  N(a),  a  g  A,  satisfies  this  assumption,  then  we  may  define  the 
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Fourier  transforms 


(3.33)  gau . . .  >  K) 

=  J  •  •  •  JeXp  |-i  £  IjUj^  rau...tttk(ux,  •  •  •  ,  Mk-i)  dux  ■  duk _ x , 

understanding  +  •  •  •  +  =  0.  For  completeness  we  set 

(3-34)  9att)  =  ra  =  qa(t)  =  pa(0, 

in  the  case  A:  =  1.  We  now  have 

Theorem  3.3.  Let  the  point  process  N(a),  a  e  A,  satisfy  Assumption  3.2.  Then 
the  process  satisfies  Assumption  2.2.  Its  cumulant  spectra  are  given  by 

(3-35)  >  ^k) 


=  (2*r‘  +  ’L  £  [FMR-ay}] 

£—  1  0-  —  l  Jev  i 

•••  cn 

jeve  jev  i 

with  the  summation  extending  over  all  partitions  (vl9  •  •  ■  ,  ve)  of  (1, 
As  particular  cases  of  the  cumulant  spectra  we  mention 


z  y. 

i^e 

k). 


(3.36) 

(3.37) 


fa  =  ra, 

f'..M  =  (2  n)~\g'aJk)  +  9a] 
=  (2  n)~^ 


J*  exp  {-ikt}  ra  a(t)  dt  +  ra 
in  agreement  with  Bartlett  [4],  p.  183.  Also 

dt 


(3.38) 

and 


fa,bU)  =  (2w)  1  J  exp  {—ikt)  ra  b(t) 


if  a  b, 


(3.39)  X2)  =  (2k)-2  [g;,„,„(A„  A2)  +  sk.tf,)  +  g;,„(A2) 

+  g'a.ai  —  ^1  —  A2)  +  g^], 

We  have  the  following  relation,  inverse  to  (3.35), 

(3-40)  9ai,-  -,ak(^l>  ‘  ‘  *  5  ^fc) 

=  t  t  -  I)-  (2lt)‘-'[n  <5{«i  -  a,}] 

^=1  <*i,  •■•,<*,.=  1  Jevi 

•  •  •  tn  <H<x,  -  Z  V  '  -  ’  Z  Ai)> 

Jev^  Jevi  jev,, 

where  the  summation  is  again  over  all  partitions  (v1}  •  •  •  ,  v^)  of  the  integers 

(1,  ••*,*). 

In  Section  2  we  discussed  a  class  linear  time  invariant  operations  on  stationary 
interval  processes.  It  may  be  of  interest  to  indicate  a  subclass  of  these  operations 
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which  carry  point  processes  over  into  point  processes.  Let  a v  j  =  0,  +  1,  •  •  •  , 
be  a  sequence  of  real  numbers.  Let 


(3.41)  a(A)  =  the  number  of  Gj  e  a, 

then. 


(3.42) 


Y(  a)  =  a*N(A ) 


=  j  a(A  —  u)  dN ( u ) 

will  be  a  real  valued  point  process  in  the  case  that  A7(a),  a  e  A,  is  one.  If  t 
j  =  0,  +  1,  •  •  •  ,  denote  the  times  of  events  of  a  realization  of  N(  a),  then  events 
of  this  Y( A)  occur  at  the  times  Tj  +  (Tk,j,k  =  0,±l,---. 

Daley  [13]  discusses  the  second  order  spectral  theory  of  point  processes,  con¬ 
siders  operations  on  point  processes,  and  presents  a  variety  of  examples. 

Example  3.3.  Suppose  that  Y(/).  —  oo  <  /  <  go,  is  an  r  vector  valued 
stationary  time  series  satisfying  the  conditions  of  Example  3.1  and  having 
Cramer  representation 


(3.43) 


Y (t)  =  J  exp  [iXt]  dZY(X). 


Suppose  A7(a),  a  €  A,  is  an  independent  stationary  point  process  satisfying 
Assumption  3.2  and  having  spectral  representation 


(3.44) 


At(a)  = 


exp  {iXt}  dt 


dZN(X). 


In  Section  6  of  the  paper  we  will  consider  the  process 


(3.45) 


X(a) 


=  1 


Y(/)  dN(t) 


=  Y(tx)  +  •  •  •  +  Y(t„), 


if  t t„  are  the  events  of  At(a)  in  the  interval  a.  One  can  check  that  this 
process  satisfies  Assumption  2.2.  If  its  spectral  representation  is 


(3.46) 


X<4)  =  J  1  exp  { iXt }  dt 


^ZX(A), 


then  we  see  directly  that 


(3.47)  dZx{X)  =  J  [dZY(X  -  a)]  dZN( a), 

for  —  co  <  X  <  oo.  Expression  (3.47)  may  be  used  to  determine  the  cumulant 
spectra  of  X(a)  in  terms  of  those  of  Y (t)  and  At(a). 

We  mention  that  Walker  suggested  the  consideration  of  real  valued  processes 
of  the  form  (3.45)  in  the  discussion  of  Bartlett  [3]. 
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4.  Stochastic  properties  of  finite  Fourier  transforms 

We  now  turn  to  an  investigation  of  certain  statistics  useful  in  the  estimation  of 
the  cumulant  spectra  of  a  stationary  interval  function  X(a).  a  e  A.  We  will  sup¬ 
pose  that  the  values  of  X(a)  are  available  for  a  contained  in  the  support  of  a 
function  h(t/T),  T  =  1,  2,  •  • •  .  We  set  down. 

Assumption  4.1.  The  function  h(t),  —  oc  <  t  <  go.  is  measurable  in  t, 
bounded,  zero  for  |^j  >  1  and  there  exists  a  finite  K  such  that 

(4.1)  J  | h(t  +  u)  —  h(t) |  c It  <  K\u\ 
for  all  real  u. 

The  inequality  (4.1)  will  be  satisfied  if  ^(/)  is  of  bounded  variation,  for  example. 
For  given  T,  the  function  h(t/T)  has  been  called  a  taper  by  Tukey  [22],  It  has 
also  been  called  a  data  window. 

The  principal  statistics  of  our  analysis  of  interval  processes  are  the  finite 
Fourier  transforms, 

(4.2)  d^\X)  =  J  ha{t/T)  exp  { -tit}  dXa(t), 

a  —  1,  •  •  •  ,  r,  —oc  <  X  <  oo.  In  the  case  of  Example  3.1,  the  statistic  (4.2)  takes 
the  form 

(4.3)  d{aT)(X )  =  J ha{t/T)  exp  {- tit }  Ya(t)  dt, 

that  is,  it  is  the  Fourier  transform  of  the  tapered  values  that  was  considered  in 
Brillinger  and  Rosenblatt  [9].  In  the  case  of  Example  3.2,  if  we  let  ifl(l),  •  *  • ,  t  a(na) 
denote  the  times  of  events  of  the  ath  sort  that  occur  in  the  support  of  ha(t/T), 
then  the  statistic  (4.2)  has  the  form 

(4-4)  J  ha(*aU)/T)e*V  {-ifoaij)}- 

j=  1 

This  statistic,  excluding  the  taper,  was  considered  in  Bartlett  [3]  for  the  case 
r  =  1  and  suggested  for  the  case  of  general  r  by  Jenkins  in  the  discussion  of  that 
paper.  In  the  case  of  Example  3.3,  the  statistic  has  the  form 

(4.5)  X  K(?j/T)exp  {-Uzj}  Ya(Tj), 

j=  i 

if  Tj,  •  •  •  ,  t„  denote  the  times  of  events  of  the  process  N( a)  in  the  support  of 
ha(t/T). 

We  next  present  a  basic  theorem  indicating  the  asymptotic  joint  cumulants  of 
the  Fourier  transform  (4.2).  In  the  theorem  we  let 


(4.6) 


Hau---,aM  =  J  Kft)  ■  ■  ■  hak(t)  exp  {-at}  dt. 
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Theorem  4.1.  Let  the  process  X(a),  a  g  A,  satisfy  Assumption  2.2.  Let  ha(t), 
a  =  \  ...  r,  — oo  <  t  <  oo,  satisfy  Assumption  4.1.  Then  as  T  — ►  go 

(4.7)  cum  ,  <P(Ak)} 

=  (2k)‘-,/„„....w(A1,  •  ■  - .  +  0(1) 

for  ax,  •  •  •  ,  ak  =  1,  •  •  •  ,  r;  k  =  1,  2,  •  •  •  .  The  0(1)  term  is  uniformly  bounded 
%n  A  | ,  *  *  *  .  Xk . 

We  see  that  the  joint  cumulants  are  of  reduced  order  unless  Xj  is  near 
zero.  We  see  from  (4.6)  and  (4.7)  that  the  joint  cumulants  based  on  disjoint 
stretches  of  data  are  of  reduced  order  as  well. 

If  ha(t)  =  1  for  0  ^  t  ^  1  and  ha(t)  =  0  otherwise,  then  this  theorem  has 
identical  nature  with  the  key  theorem  used  in  Brillinger  and  Rosenblatt  [9], 
Brillinger  [6],  Brillinger  [7]  to  develop  properties  of  spectral  estimates.  The 
results  of  these  papers  therefore  become  directly  available.  We  indicate  a  selection 
of  results  that  now  hold. 

We  begin  by  considering  the  asymptotic  distribution  of  the  finite  Fourier 
transform.  Let  Nf(p,  X)  denote  the  complex  r  variate  normal  distribution  with 
mean  /i  and  covariance  matrix  £.  We  have 

Theorem  4.2.  Let  X(a),  a  g  A,  be  an  r  vector  valued  interval  process  satisfying 
Assumption  2.2.  Let  Sj(T)  be  an  integer  with  Xj(T)  =  2itSj{T)/T  -*■  Xj  as  T  -*■  oo 
for  j  —  1,  •  •  •  ,  J.  Suppose  Xj(T)  ±  Xk{T)  0  for  j,  k  =  1,  ■  ■  •  ,  J.  Let 

(4.8)  d^T)(A)  =  f  exp  {  —  iXt]  dX{t) 

Jo 

for  -oo  <  A  <  oc.  Then  d (xT){Xj(T)),  are  asymptotically  indepen¬ 

dent  iV^(0,  2nTi'x  x{X))  variates,  respectively .  Also  dxT)(0)  =  X(0,  71]  is  asymp¬ 
totically  iVr(7T^,  2nTf'x ,x(0))  independently  of  the  previous  variates. 

This  theorem  has  the  nature  of  a  central  limit  theorem.  Let  Wf{n,  £)  denote 
the  complex  Wishart  distribution  of  dimensions  r  x  r,  degrees  of  freedom  n 
and  covariance  matrix  X.  Define  the  matrix  of  periodograms 

(4.9)  r/W)  =  (2*7’)- 1  d‘/>(A)  Wvi'- 

We  have  the  following  corollary. 

Corollary  4.1.  Under  the  conditions  of  Theorem  4.2,  if  Ax  =  •  •  •  =  Xj  =  A 
and  if 

(4.10)  =  t  UMUO). 

i=  1 

fx'fx  (A)  is  asymptotically  J~  1  Wf{J.  f’x,  x  (A))  as  T  — >  oo . 

This  corollary  makes  precise  the  chi  square  approximation  for  the  distribution 
of  second  order  spectral  densities  of  point  processes  suggested  by  Bartlett  [3]. 
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We  next  construct  consistent  asymptotically  normal  estimates  of  the  cumulant 
spectra  of  different  orders  of  an  interval  process  X(a),  a  €  A.  We  begin  by  letting 
W(ux,  •  •  •  ,  uk)  be  a  weight  function  satisfying 

Assumption  4.2.  The  function  W{ux,  •  •  •  ,  uk ),  —  x  <  Uj  <  x,  is  symmetric 
in  ux,  •  •  •  ,  uk,  is  concentrated  on  the  plane  X*  Uj  =  0,  and  is  such  that 


(4.11) 


W I 


<5(  Uj  dux  ■  •  ■  duk  —  1 


and 


(4.12) 


W  (  tlj,  •  ,  Uk_x, 


■  Uk-  1  • 


g  A 


I 


i 


-k -£+  1 


for  some  A,  £  >  0,  /  =  1,  •  •  •  ,  k. 

Given  the  sequence  of  nonnegative  numbers  BjK  T  —  2,  3,  •  •  •  .  we  set 
(4.13)  WT(u1,- ■  ■  ,uk)  =  (Bif)~k  +  1  W(B{f)"u1,  •  •  •  ,B{fr'uk). 


We  suppose  B(f]  ^  B(f]  ^  .  Next  we  set  •  •  •  ,  uk)  =  1  if  L*  Uj  =  0 

but  no  proper  subset  of  the  Uj  has  sum  0.  and  set  it  =  0  otherwise.  Let 

(4.14)  d(xr)(A)  =f  exp  {  —  ikt)  dX{t). 

Jo 

Finally  set 

(4.15)  K)  =  (2n)-»*,7’-'  n  <>(A j). 

j=  1 

As  an  estimate  of  fai  . . .  >0k(Ai,  •  •  •  ,  Ak)  we  now  take 
(4-16) 


■  I  H  V 


2ns  j 

~Y~ 


k)c- 


In  connection  with  this  estimate  we  have  the  theorem, 

Theorem  4.3.  Let  X(a),  a  e  A,  satisfy  Assumption  2.2.  Let  W{ux,  •  •  •  ,  uk) 
satisfy  Assumption  4.2.  Let  f{aTf ...  ak(A1,  •  •  •  .  kk)  be  given  by  (4.16).  Let  B{£]  —>  0, 
W)‘  lT—>ooasT—>co,  then 

(4.i7)  fi/:r.,.....,(A,.---,Al) 
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=  J  •  •  •  J  HWi  -  a,,  •  •  •  ,  Xk  -  a •  •  •  -  a*) 

- ^(a,  4-  •  •  •  +  afc)  dcti  ■  ■  ■  dctk  +  0(B§y~lT) 

=  4)  +  o(B»>)  +  o(B«»"r), 

(4.18)  to  (fiWf-T  Cov  ,  4). •  • ,  ft)} 

=  -71  X  ‘H^l  ~  Pp,  i}  —  fiP,k}faiap ,  1  (^l)  "  f akaP'  k  (^k) 


where  the  summation  is  over  all  permutations  P  of  the  integers  l,  ■■■,  Ic.  Collections 
of  spectral  estimates  are  asymptotically  jointly  normally  distributed  as  T  — ►  oo 
estimates  of  different  orders  asymptotically  independent  and  estimates  of  the 
same  order  having  covariance  structure  given  by  (4.18). 

We  next  turn  to  the  development  of  an  empirical  analysis  of  the  linear  time 
invariant  model, 

(4.19)  Y(a)  =  a*X(A)  +  8(a) 


=  J  a(A  —  u)  dX(u)  +  e(A), 


with  X(a),  8(a),  a  e  A,  independent  stationary  interval  processes  and 


J  M  tf|a(«)| 


In  differential  notation  we  may  write  (4.19)  as 


dY(t)  =  J  \d&{t  —  «)]  dX(u)  +  de(t). 


Denote  the  cross  spectral  density  matrix  of  the  process 


~X(  A  )~ 
_Y(a)_  ’ 


a  €  A,  by 


fx.xU)  fx.rUr 


and  that  of  8(a),  a  e  A.  by  f'E  E(A).  Set 


A(;')  -  J  exp  {  —  ikt]  da,(t). 
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Then  (4.19)  gives 

(4.25)  fy,v(^)  =  A(A)fx>  X(A), 

(4.26)  fy>r(A)  =  A(A)f^sX(A)A(I)t  +  fe.,(A). 

These  last  suggest  that  we  may  base  estimates  of  A(k)  and  f'  e(A)  on  an  estimate 
of  the  spectral  density  matrix  (4.23).  We  could  construct  an  estimate  of  this 
last  in  the  manner  of  (4.16);  however,  in  order  to  display  an  alternate  form  of 
spectral  estimate  of  order  two  wTe  proceed  slightly  differently. 

In  constructing  this  alternate  estimate  we  let  h{t),  —  oo  <  t  <  oc .  be  a  tapering 
function  satisfying  Assumption  4.1.  We  then  set 

(4.27)  //fc(A)  —  j  h(t)k  exp  {  —  ikt]  dt. 
for  —  oo  <  k  <  oc.  We  next  set 


d(/J(A)  =  J  h(t/T)  exp  {  —  ikt)  dX(t). 

(4.28) 

d(/*(A)  =  I  h(t/T)  exp  {  —  ikt)  dY(t) : 

and  we  let  W(ct)  be  a  weight  function  satisfying 

Assumption  4.3.  IF(a),  —  oc  <  a  <  oo  .is  real  valued ,  even,  absolutely  inte¬ 
grate.  has  an  absolutely  integrable  first  derivative,  and 

(4.29)  P  W(<x)doi  =  1. 

J  -  00 

The  variate  (4.22)  has  mean 

(4.30)  l*  j4|. 

Y  J 

Estimates  of  f*,  f'K  based  on  tapered  values  are  provided  by 

fU  =  J  h(t/T)dX(t)/j  h(t/T)  dt  =  dY(0)/[ra,(0)], 

(4.31) 

f(,T>  =  j  h(t/T)  dY(l)/j  h(t/T)dt  =  dVr>(0)/[7’//,(0)], 

respectively.  The  Fourier  transform  of  the  process  (4.22)  corrected  for  its  sample 
mean  is  then  given  by 

e(/)(; l)  =  J  exp  {-ikt)h{t/T)[dX(t)  -  dt] 

=  d</>(A)  -  d(P (0)H Y{T k)/ H j(0), 


(4.32) 
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ei,T'(A)  -  J  exp  {-iXt}h(.t/T)[dY(t)  -  f}TI  d(\ 

=  d'/’U)  -  4T|(0)//1(7'A)///1(0). 

Let  Bt  be  a  sequence  of  nonnegative  numbers  tending  to  0  as  T  -*■  oo .  Set 
(4.33)  W(T){<x)  =  Bj 1  W(Bj  1ol). 


As  an  estimate  of  the  cross  spectral  density  matrix  (4.23)  we  now  propose 


(4.34) 


fx\  x  (^) 

fyfx  (^) 


fxT r  (^)  _  f 

fy[y(^)_  ^ 


Wm{X  -  ct)  (InTy 1 


As  estimates  of  A(A),  f'£(A),  we  then  take 

(4.35)  A(T)(A)  =  f  aa)fB(A)-1, 

(4.36)  g%(A)  =  fRvw  -  f^(A)f 

We  can  now  state  the  following  theorem. 

Theorem  4.4.  Let  the  process  X(a),  a  g  A,  satisfy  Assumption  2.2.  Suppose 
fx  x(^)  nonsingular.  Let  the  process  «(a),  a  g  A,  satisfy  Assumption  2.2.,  have 
mean  0  and  he  statistically  independent  of  the  process  X(a),  a  g  A.  Let  a(A)  satisfy 
(4.20).  Let  Y(a)  be  given  by  (4.19).  Let  W(cc)  satisfy  Assumption  4.3  and  h{t) 
satisfy  Assumption  4.1.  Then  if  BT  — ►  0,  BtT  -*■  co  as  T  -*■  go, 

(4.37)  lim  a"ve  A(T)(A)  =  A(A), 


(4.38)  lim  BtT  Cov  {vec  A(T)(A),  vec  A{T)(p)} 

T  -*  oo 

=  2jtH4(0)ff2(0)-2d{A  -  /j}r,.,(A)  ®.f'x,xar'  J  W(a)2da, 

(4.39)  lim  ave  g 2(A)  =  f',c(A), 

T  ->  oo 

(4.40)  lim  BTr  Cc&  {fl$(A),  gU?„Ui)} 

r  -►  oo 

=  2nH4(o)fi2(or2(d{A  - 


(4.41)  lim  BtT  Cov  {vec  A(T)(A),  g^lifi)}  =  0, 

T~*  oo 

for  j,  k,  m,  n  =  1,  ■  "  ,  *•  */&e  variates  A(7)(A),  ^(/i)  are  asymptotically 

jointly  normal  with  the  above  covariance  structure. 
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(In  this  theorem  ave,  Cov  have  technical  definitions  allowing  the  use  of  Taylor 
series  expansions  in  determining  asymptotic  moments.  See  Brillinger  and  Tukey 
[10].) 

In  the  case  r  =  s  =  1,  we  may  define  the  coherency  of  X(a)  with  Y( a)  at 
frequency  X  by 


(4.42) 


^Y,x(^)\2 


Ui.xWi.rM]' 


As  an  estimate  of  the  coherency  we  consider  the  statistic 


(4.43) 


KW)  I2 


l/«U)|2 

UkW)fiW)T 


We  then  have  from  the  theorem 

Corollary  4.2.  Under  the  conditions  of  the  theorem  and  iffx.x  (/0  /y,  y  (^)  ^  0, 
|-^yT)v(^-)|2  is  asymptotically  normal  with 


(4.44)  lim  a"ve  |.fty7\(A)|2  =  \Ry  x(^)\2 

r->  oo 

and 


(4.45)  lim  BtT  Cov  { | (A)|2.  K’xMl2} 

T-»  oo 

=  4nH 4(Q)H 2(0)~ 2  \_S{X  -  p] 

+  <5{A  +  /i}]  |/?y  x(A)|2[l  —  | Ry, x (^-) | 2]2  j  H(a)2  da. 

A  comparison  of  the  results  of  this  theorem  and  its  corollary,  with  the  corre¬ 
sponding  results  for  the  regression  of  one  vector  valued  stationary  time  series 
on  another,  shows  that  they  are  identical.  This  will  also  be  the  case  for  the 
interval  process  extension  of  many  of  the  asymptotic  results  of  the  analysis  of 
stationary  time  series. 


5.  Estimation  of  product  densities 

Let  N(a),  a  g  A.  be  a  stationary  point  process  satisfying  Assumption  3.2.  We 
have  defined  various  characteristics  of  such  a  process.  These  may  be  summarized 
as  follows : 

(5-1)  ?ai,..,flk(li,"',4) 

=  lim  p{dN  {tU  =  1,  •  •  •  ,  dN  (tk)  =  l}/(dt1  •  •  ■  dtk) 

for  tx .  •  •  •  .  tk  distinct ; 

(5.2) 

=  Z  -  U'PayJevftjJ  e  V1  )  '  *  '  PayJ^Uj' J  G 

r=i 


INTERVAL  FUNCTIONS 


503 


(5.3) 


(5.4) 

(5.5) 


Pai,  ■  •  ■ ,  ak  Ol  ’  ’  ^k) 

k 

=  Z  E  Vl)  '  '  '  E 

l 


rai,---,«k(“l»  ’  •  ’  »  Mfc-l)  =  0«,.---.ak(“l>  ’  '  ’  -  Uk-  1-  0)  ; 

9ai,  •  ■  ■  .ajcC^l  •  ’  ^Ic—  l) 


)  diii  •  •  ■  duk_1 ; 


and  if  ,  Ak_j)  is  integrable 


(5.6) 

•  •  •  ,uk- j) 


=  (2w)  fc+1  J'"  J  exP  j*  Z  “  dK-i- 


Also 


(5-7)  9au  -,ak(^  1>  >  ^k)  ~  9ai,-  ■  ,ak(^l  >  ^k-l) 

understanding  £*  =  0.  Continuing 

(5-8)  fai,...,ak(ln--- ,h) 

=  (2n)-k+i  z  z  cn  ^(«i  -  aii]  •  •  •  cn  h**  -  °ji] 

^=1  ai,-..,a^=l  Jev  i  Jev, 

•^„  ..,a,[Z  4  ‘  •  Z  ;J: 

Jev !  jev  i 

(5.9)  9ai,  •  ■  ■ ,  ak  (^1  >  >  ^ k ) 


=  z  z  (-irv-i)![n<H«i  -  «;■}]•• 

<?=1  a  i ,  ■  •  • ,  a,  =  1  Jevi  Jev, 

I  hi 

Jev  i  Jev, 


The  summations  in  (5.2),  (5.3),  (5.8),  (5.9)  are  over  all  partitions  (Vj,  •  •  •  ,  v,,)  of 
the  integers  (1,  •  •  •  ,  k). 

In  the  previous  section  we  developed  an  estimate  of /ai,...,0k(^i>  ‘  '  *  »  K)-  Let 
us  now  put  this  work  to  use  in  developing  estimates  of  the  various  functions 
listed  above.  As  an  estimate  of  gau. . .  ,ak{^ i-  ' '  *  ,  Afc).  in  the  light  of  (5.9),  we 
may  consider 

(5.10)  fltf?.... 

=  i  t  (-ii'-v-i).[n«{«1-./}]-[n^-(i] 

(  =  1  ai,  •  •  ■  ,  a,  —  1  Jevi  jevr 

•12-ti*  r,r  ,ii  ',■■■■■!.  hi 

jev  i  jev , 
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where  the  (pl5  •  •  •  ,  fie).  /  =  1.  •  •  •  ,  k,  are  formed  in  the  manner  of 

Theorem  4.3.  From  that  theorem  we  see  that 

(5.11)  ,h) 

=  +  0(/*?>)  +  OW'r1); 

and  because  estimates  of  order  less  than  k  have  asymptotic  variance  of  smaller 
order  than  that  of  estimates  of  order  k,  the  covariance  of  g^].  ■ .  <ak(k1.  '  *  '  •  Ak) 
with  .  ,!,(/*  i,  •  •  •  ,  fa)  will  be  asymptotically  equivalent  to  that  of  y(tk 
(A1?  •  ■  •  ,  kk)  with  bf  (/i  i,  •  •  •  ,  Ht)  as  given  in  Theorem  4.3.  Also  the  estimates 

will  be  asymptotically  normal  and  estimates  of  different  orders  will  be  asymp¬ 
totically  independent. 

Suppose  next  that  g'ai  ... •  •  •  ,  Ak_j)  vanishes  for  \kj\  >  A.  As  an  esti¬ 
mate  of  rai  ...  ak(u i,  •  •  •  ,  wk_  i)  we  can  then  consider 

(5.12) 

=  (2tc)"*+1^— Jexpj*  X  -  X  ^ 

d  k  ^  _  j . 

From  (5.11)  this  estimate  will  be  asymptotically  unbiased. 

By  analogy  with  Theorem  5.2  of  Brillinger  [6],  we  would  expect,  for  example, 
that 

(5.13)  lim  T  Cov{r(Z]bx{ux),r(J2]b2{u2)} 

1  —*  cc 

=  J^AeXP  “  U2  )}faua2((X)fbl.b2(-^)  dci 

+  J_AeXp  ^a(Ml  +  M2)}/o!,b2(a)/bI,a2(-a) 

+  2n  jj:  exp  {i(a1u1  +  a2M2)}/;it6l>a2>b2(a1,  -oc^.  a2)  dotlda2: 


in  the  case  k  =  2  and  aj  =/=  bj. 
Next  one  can  take 


(5.14)  qZ\. =  *•£.’.. .,„(<!  -  -«*) 

as  an  estimate  of  qau . . .  <ak{t 1?  •  •  •  ,  tk)  and 


(5.15) 


pZI 


Mi,  •••.«*)  =  X  dfi&Mj'J  E  vi)  •  •  ‘  e 


e  =  i 


as  an  estimate  of pai  ... >flk(<1,  •  •  •  ,  fk). 

In  the  case  k  =  1,  we  would  estimate  ra  by  iVa(0,  T~\jT .  In  Theorem  4.2  we 
saw  that  this  statistic  was  asymptotically  normal  with  mean  ra  and  variance 

2xT-%.(0). 


INTERVAL  FUNCTIONS 


505 


6.  Estimation  of  second  order  spectra  from  sampled  values 

Let  Y(t),  — oo  <  t  <  oo,  be  a  real  valued  time  series  satisfying  the  conditions 
of  Example  3.1,  having  mean  c'y  and  autoeovariance  function  cYY(u),  —  oo 
<  u  <  oo.  Let  N( a),  a  g  A,  be  an  independent  real  valued  point  process  satis¬ 
fying  Assumption  3.2,  having  mean  intensity  rN  and  autoeovariance  density 
rNN(u),  —  oo  <  u  <  oo.  Suppose  that  events  of  a  realization  of  the  process  N(  a) 
occur  at  the  times  ■  •  •  ,  r„  in  the  interval  (0,  T~\.  Consider  the  problem  of  esti¬ 
mating  the  autoeovariance  c'y  y(w),  —  oo  <  u  <  oo,  and  power  spectrum/'^  y(A), 
—  oo  <  A  <  oo,  of  the  series  Y(t)  from  the  values 

(6.1) 

and 


(6.2)  ¥(!,),■  ■•  ,Y(zn). 

We  can  construct  a  stationary  interval  process  X(a),  a  e  A,  in  the  manner  of 
Example  3.3  by  setting 

(6.3)  X(A)  =  f  F(l)^(l), 

or,  in  differential  notation,  by  setting 

(6.4)  dX(t)  =  Y(t)dN{t). 

The  first  and  second  order  measures  of  this  process  satisfy 

(6.5)  Cx  dt  =  c'YrN  dt, 
and 

(6.6)  dC'x  x(u)  dt  =  (c'YiY{u)rNiN(u)  +  c'XY(u)rNS(u)  +  c'y>y(M)r£ 

+  ( c'Y)2rN,N(u )  +  (c'Y)2rNS(u))  du  dt. 

The  measure  C'x <x(u)  is  seen  to  have  absolutely  continuous  part  and  an  atom 
of  mass  c'y  y(0)rN  +  ( c'Y)2rN  at  u  =  0.  If  we  let  rx <x(u)  denote  the  derivative  of 
the  absolutely  continuous  part  of  C'x<x(u)  then,  from  (6.6), 


(6.7)  rx,x(u)  —  cy,  r (u)rN, n(m)  4"  cy  y(w)r^  +  (cy)  rN,N(u) 
for  —  oo  <  u  <  oo .  For  convenience  set 

(6.8)  h(u)  =  rx,x(u )  ~  (c'Y)2rN,N(u)- 

If 

(6.9)  +  rN,N(u )  ^  6, 
then,  from  (6.7), 


h{u) 

[rN,/v(M)  +  rN~] 


(6.10) 


cy,y(u)  — 
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We  see,  from  (6.8)  and  (6.10),  that  an  estimate  of  c'Y  Y{u)  may  be  constructed 
from  estimates  of  rxx(u),  c'Y,  rN  N(u ),  rN.  One  can  then  proceed  to  form  an  esti¬ 
mate  Of /y  y(A). 

Alternatively  we  could  proceed  directly  to  the  frequency  domain  and  note 
that  the  power  spectrum  of  the  process  A(a)  is  given  by 


(6.11)  fXJS(k)  =  (2k)-1 


rx  x(u)  exP  { — du  +  c'yj(0)rN  + 


—  a)9rv,v(a)^a  +  /y,y(^-)rv  +  (cy)2/jv,  jv(^) 

+  (271)  1Cyy(0)%, 

for  —  x  <  A  <  x  .  If  we  rewrite  this  in  the  form 

(6.12)  /y,y(^-)^  +  j  fy,  y (a)9,N,  —  &■)  da 

=  f x,x(^)  ~  (cy)V n,n(^)  ~  (2?r)  1Cy  y(0)rN 

=  i^(A), 

then  we  have  an  integral  equation  for/y  y(/l).  This  equation  may  be  solved  for 
/y  y(A),  under  the  condition  (6.9),  as  follows:  set 

(6.13)  P(X)  =  (27t)_1  J  exp  {  —  iXu}rNN{u)/\_rj,  +  rN  N(u)~\  du, 
then 

(6.14)  /y  y(A)  =  r^2H(X)  -  2nrx2  J  ,P(A  —  a)H(a)da. 

Once  estimates  of  rN  N(u),  c'Y.  c'Y  Y{0),  f'N  N{X).  fx  x{X),  are  available  an  esti¬ 
mate  of/y  y(A)  may  be  constructed  from,  (6.14).  The  estimates  may  be  determined 
as  follows : 

(6.15)  r(P  =  n/T : 

(6.16)  4T)  =  [7(10  +  •••  +  Y(T„)]/n: 

(6.17)  <V(0)  =  [7(t,)2  +  •  •  •  +  Y(T„)2]/n ; 

(6.18)  c(yrV(0)  =  m(yrV(0)  -  c(yr)2; 

and  finally  estimates /^(^b/^xU)  may  he  constructed  in  the  manner  of  (4.16) 
or  (4.34). 

A  problem  related  to  the  one  just  considered  is  that  of  obtaining  as  estimate 
of  the  cross  spectrum  fY  t  y2(  X)  of  a  series  Yi  (t)  with  a  series  Y  2{t)  from  the  values 

b  l(Tlb  ■ .  J  b  i(t„), 


(6.19) 

(6.20) 
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and 


(6-21)  Y2(x1)r-tT2(xa). 

In  this  case  the  expression  (6.11)  is  replaced  by 


(6.22)  fxi,x2(^)  —  j  f,Y1,Y2( ^  ~  a)9rN,iv(a)  +  f'YuY2(^)rN 

+  (CYi)(CY2)f  'N,  n(^)  +  (2JC)  1Cyljy2(0)rJv. 

A  second  related  problem  would  be  to  construct  an  estimate  oif'Ylty2(X)  from 
the  values 


(6.23) 

(6.24) 

(6.25) 

(6.26) 


Gl )  "  '  i  ffmi 

Tl,  *  *  *  ,  T„, 

Fi((Ti),  •  ,  Yx(om), 
Y2(t: i),  '  ‘  ,  Y2(zn), 


where  ol,  •  •  •  ,  <rm  are  the  times  of  events  in  (0,  T]  of  a  point  process  N a)  and 
Tj,  •  •  •  ,  t„  are  the  times  of  events  in  (0,  T~\  of  a  related  point  process  N  2{  a)  with 
the  bivariate  point  process  satisfying  Assumption  3.2.  In  this  case  expression, 
(6.11)  is  replaced  by  the  simpler  expression 

(6.27)  f  'xi,x2(^)  =  J  Iyi.yA^  ~  <x)9Ni,N2(<*)  +  /yl(y2(^)rJV,rJV2 

+  cYicyJnun2{ ^)- 


7.  Further  considerations 

We  next  discuss  briefly  some  practical  implications  and  extensions  of  the 
previous  results.  We  saw,  in  Section  2,  that  if  X(a),  a  e  A,  was  a  stationary 
interval  process  with  cumulant  spectra 

(7-1)  fxai,-,xaUi,---  ,h), 

then 

(7.2)  Ya(t)  =  J  <f>a(t  -  u)  dXju), 

a  =  1 ,  •  •  •  ,  r,  —  oo  <  t  <  oo,  was  a  stationary  time  series  with  cumulant  spectra 

(7.3)  /y0l,-",rak(^i’  ‘  *  *  >  h)  =  ^><i1(^i)  '  ’  *  ®ak(^k)fxai>  - ,xak(^ i»  '  *  ‘ 

This  suggests  that  one  might  estimate  the  spectrum  (7.1)  by  a  statistic  of  the 
form 


(7.4) 
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having  formed 

(7-5)  y.tf„  ••••■**) 

in  the  manner  of  Brillinger  and  Rosenblatt  [9].  (in  the  case  k  =  2  this  sug¬ 
gestion  was  made  by  Priestly  in  the  discussion  of  Bartlett  [3].)  This  procedure 
is  seen  to  be  analogous  with  the  technique  of  prewhitening  a  time  series  prior 
to  estimating  its  spectrum.  This  analogy  suggests  that  we  should  choose  the 
(fta(X)  so  that  the  spectrum  (7.3)  is  near  constant  for  Xj  in  some  finite  region. 
The  estimate  (7.4)  is  seen  to  have  the  important  advantage  of  allowing  the  use 
of  existing  spectral  programs  and  also  of  allowing  a  simultaneous  prewhitening 
of  the  data. 

The  proposed  analysis  may  be  related  to  the  analysis  of  a  continuous  time 
series  in  another  way.  The  basic  statistic  of  our  analysis  is 

(7.6)  f  exp  {  —  iXt}h{t/T)  dX(t). 

Jo 

If  we  approximate  (7.6)  by  a  Stieltjes  sum,  then  we  obtain 
r-  i 

(7.7)  Y  exP  {  —  iXt}h(t/T)[X(t  +  1)  —  X(/)]. 

t  =  o 

An  examination  of  expression  (7.7)  shows  that  it  corresponds  to  carrying  out  an 
empirical  spectral  analysis  on  the  time  series  of  first  differences.  This  procedure 
is  common  in  the  analysis  of  economic  time  series. 

Computations  involved  in  forming  (7.6)  may  be  prohibitive.  Therefore  there 
is  much  to  be  said  for  a  procedure  involving  splitting  the  data  into  AT  segments 
of  length  S,  forming  an  estimate 

(7.8)  /£’... 

for  the  wth  segment,  n  =  1,  •  •  •  ,  N ,  and  taking 

(7-9)  A"1  t  f If!  aMf-'--K)n 

n=  1  ’ 

as  a  final  estimate.  Authors  recommending  such  a  procedure  include :  Bartlett 
[2],  Welch  [23],  Lewis  [19],  and  Huber  et  al  [15].  The  asymptotics  of  such 
estimates  are  directly  determinable  from  the  results  of  Theorem  4.3  because, 
following  the  remark  after  Theorem  4.1,  Fourier  transforms  based  on  disjoint 
stretches  of  data  are  asymptotically  independent.  A  variety  of  further  remarks 
concerning  practical  aspects  of  the  calculations  in  the  case  of  a  point  process  are 
made  in  Lewis  [19]. 

We  remark  that  the  calculations  proposed  in  this  paper  reduce,  in  the  case  that 
the  interval  process  X(a),  a  e  A,  is  an  integral  of  a  continuous  time  series,  to 
the  usual  calculations  of  the  frequency  analysis  of  time  series. 

Extensions  of  the  definitions  and  theorems  of  this  paper  to  a  case  in  which  t 
is  vector  valued,  t  e  Rp,  appear  fairly  immediate  if  one  takes  the  approach  of 
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Brillinger  [7].  A  different  sort  of  extension  would  result  from  a  consideration  of 
processes  whose  differences  of  higher  order  than  the  first  are  stationary  (see 
Yaglom  [24]). 

8.  Proofs 

Proof  of  Lemma  2.1.  If  Mai  ...  ak(tx,  ■  ■  ■  ,.tk)  corresponds  to  the  measure 
determined  by  the  coordinates  tx,  •  •  •  ,  tk,  let  Nai  ... t0k(ux,  •  •  •  ,  uk_x,  tk)  corre¬ 
spond  to  the  measure  determined  by  the  coordinates  ux  =  tx  —  tk,  •  •  •  ,  uk_  x  = 
tk_  j  —  tk,  tk.  The  initial  measure  is  invariant  under  the  transformation 
t  ,  •  •  •  ,  tk  ->  tx  +  t,  •  ■  ■  ,  tk  +  t.  The  second  measure  is  therefore  invariant 
under  the  transformation  tk  — ►  tk  +  t.  We  see  therefore  that 

(8.1)  Nau...,aM  1’  "  *  »  Uk-\ .  tk)  -  Nai,---,ak(u  1,  •  •  •  >  Uk_x,  0) 

=  Nai,---,ak(u i,  *  '  *  »  k  +  t)  -  Nau...tak(ux,  ■  •  •  .  Uk_x,  t). 

Suppressing  ax,  ■  ■  ■  ,  ak,  ux,  ■  •  ■  ,  uk_  x  this  last  may  be  written 

(8.2)  N(tk  +  t)  =  N(tk)  +  N(t)  -  N (0). 

Under  the  given  conditions,  all  solutions  of  this  functional  equation  have  the 
form 

(8.3)  N(tk)  =  M'tk  +  AT(0), 
giving  the  indicated  result. 

Proof  of  Theorem  2.1.  Assume  the  results  of  Theorem  4.1  hold.  It  will  be 
proved  later.  Set 

(8.4)  d(xT)(/l)  =  J  exp  {—  iXt)  dX(t), 

using  the  notation  of  Section  4  with  h{t)  =  1  for|f[  ^  1  and  h(t)  =  0  otherwise. 
One  has  therefore 

(8.5)  Z(XT)(A)  =  (27i)_  1  f  d{/]{a)da. 

J  o 

One  now  uses  expression  (4.7)  to  see  that 

(8.6)  E\Z^(X)  -  Zf>(A) \2k  -+  0 

as  8,  T  -»  oo  for  k  =  1,2,-**  ;  a  =  1,  •  •  •  ,  r.  It  follows  that  there  exists  Zx(X) 
such  that  Z(/)(A)  ->  ZX(A)  in  mean  of  order  v  for  any  v  >  0. 

One  next  checks  that 


as  T  — *  oo,  again  using  expression  (4.7).  This  gives  (2.27). 
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Finally  one  checks  that 


(8.8) 


E \Xa(A) 


exp  {ikt]  dt 

A 


dZ'J'(X)  |2 


0 


as  T  -*■  oo.  This  gives  (2.28). 

Proof  of  Theorem  3.1.  We  first  state  and  prove  a  lemma. 
Lemma  8.1.  With  the  conditions  and  notation  of  Assumption  3.1, 


(8.9) 


lim  I A  i 

lA.I-0  1 


M  (4| )]  •  •  • 

[A’« _ ( A )  •  ■  •  A-„m>  ( At )]  j 


r 

-  i 


n  siai  -  “ii 

j=  i 


mk 


n 

S{aj  -  ak} 

P*,.--.ak(*l-  •  ’  •  ,  4) 

_  j  =  mk  -  i  + 

i 

uniformly  in  tx,  •  •  •  ,  tk 

for  integers  1  ^ 

m  j  <  m2  < 

<  Wfc-i  <  mk- 

Proof 

Suppose  first  that 

■  =  «!• 

(8.10) 

amk 

-  ,  +  1  •  '  ‘ 

Atmk  =  V 

Now 

(8.11) 

-mk  - 

} 

=  I  •  • 

flj  >  J 

« 

i _ i 

Ai)  =  n,.  •  • 

'  •  ^ak(A*)  = 

=  ^[^.(Ai) 

=  1,  • 

with  the  second  summation  extending  over  some  n}  ^ 

2  and  with  | L{nx.  •••.%)!  fS 

Kd\ aa  |"1  •  •  •  |  Ak|"k  from  (3.8),  and  so 

(8.12)  :  Hmo  |a,|-‘  ■  •  •  |4t|-1£{ff„l(A,)-  ■  ■  • 


Poll .  •  •  •  »  #k  (  ^  1  ■  ’  ^k)i 

uniformly  in  tx,  •  •  •  ,  tk  from  (3.9).  Continuing  if  (8.10)  is  not  satisfied  for  some 
al,  •  •  •  ,  ak,  then  one  can  see  from  (3.8)  that  the  limit  in  (8.9)  is  0  uniformly  in 
tlf  •  •  •  ,  tk.  This  completes  the  proof  of  the  lemma. 

Turning  to  the  proof  of  the  theorem ;  let  <pj(t)  be  continuous  is  some  interval 
of  A  and  0  elsewhere  for^'  =  1,  •  •  •  ,  k.  We  have 
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By  bounded  convergence, 


(8.14)  fijj  <*>,(()  dNai(t)  ■  •  •  J  M)  dNak(t) 


=  lim  £  •  •  •  £  4MM)  •  •  •  ite  +  e]  •  •  •  Nak{iks,  ike  +  e]} 

£-0  i,  ik  K 

=  iimz  z  tn  siaj  -  ai)]  •••cri  Haj  -  a<r)] 

8-^0  ^=1  j],  —  ,  a  =1  jevi  jev,f 

•z  ■  ■  ■  z  cn  ^(*i£)]  •  •  •  cn  tAWip*  •  •  • , 

*i  if  7'evi  jev^ 

where  the  summation  extends  over  partitions  (v1?  •  •  •  ,  v,)  of  (1,  •  •  •  ,  k)  if  we 
separate  out  terms  in  (8.14)  with  the  same  argument  and  use  Lemma  8.1.  We 
now  see  that  expression  (8.14)  equals 

(8-15)  Z  Z  til  Siaj  ~  ai}]  •  *  *  HI  Siaj  ~  a)l 

\  a  !,•••,  0^  =  1  jevt  je\e 

■  f "  ■  fen  i)]  ■  ■  *  cn  •  •  •  ^t)dxx  •  •  -dit. 

J  •>  jev  i  je\f 

Expression  (3.13)  now  follows  from  (8.15)  taking  the  (f>j{t)  to  be  indicator 
functions. 

Proof  of  Theorem  3.2.  One  proves  (3.28)  from  (8.15)  and  then  obtains 
(3.22)  by  taking  the  to  be  indicator  functions. 

Proof  of  Theorem  3.3.  This  follows  directly  from  (3.28). 

Proof  of  The5rem  4.1.  Let  A*,r)(f)  =  ha(t/T).  The  cumulant  at  issue  is 
given  by 


(8.16)  J--- 


‘  exp  |  — *  X  Ajtjj  dC'au...  ak(tx  -  tk,-",tk- 1  —  h)  dtk 

■i-fii.  h<aTl)(u1  +<)••'  KV-  ,(«*-!  +  OK ?(0 


exp  <  —  i  Z  A.jt  >  dt 


k-  1 


exp  ^  —  i  £  kjUj\dC'au ...,ak{uu  •  ••  ,  «*_*). 
i 


The  indicated  result  now  follows  as 
(8.17) 


£  MVi  +  0  •  ’  • 

~  h™{t)  •  ■  •  h™{t)  exp  { -  ikt }  dt 
for  some  finite  C  following  Assumption  4.1. 


^  c  Z  bl> 
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Proof  of  Theorem  4.2.  This  follows  directly  from  Theorem  4.1  in  the 
manner  of  corresponding  results  in  Brillinger  [7]  and  Brillinger  [8]. 

Proof  of  Theorem  4.3.  This  follows  directly  from  Theorem  4.1  in  the 
manner  of  the  principal  theorems  in  Brillinger  and  Rosenblatt  [9] . 

Proof  of  Theorem  4.4.  This  follows  directly  from  Theorem  4.1  in  the 
manner  of  corresponding  results  in  Brillinger  [7]  and  Brillinger  [8]. 
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